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Abstract. The Barrett-Crane models of Lorentzian quantum gravity are a family of spin foam 
models based on the Lorentz group. We show that for various choices of edge and face amplitudes, 
If} • including the Perez-Rovelli normalization, the amplitude for every triangulated closed 4-manifold is a 

■ non-negative real number. Roughly speaking, this means that if one sums over triangulations, there 
' is no interference between the different triangulations. We prove non- negativity by transforming 

^v^j i the model into a "dual variables" formulation in which the amplitude for a given triangulation is 

expressed as an integral over three copies of hyperbolic space for each tetrahedron. Then we prove 
that, expressed in this way, the integrand is non-negative. In addition to implying that the amplitude 
is non-negative, the non-negativity of the integrand is highly significant from the point of view of 
numerical computations, as it allows statistical methods such as the Metropolis algorithm to be used 
, for efficient computation of expectation values of observables. 

<n : 

7— I . 

1. Introduction 

o ; 

Spin foam models provide a background independent approach to quantum field theory and arise 

■ naturally as a path integral formulation of quantum gravity ^ HI El 03 113 120] • 
This paper deals with the Lorentzian Barrett-Crane model, but to set the stage we first give some 

background on the Riemannian model. The Riemannian model is not expected to be physically real- 
istic, but has been an important arena for testing theoretical and computational aspects of the spin 
foam program. Shortly after Barrett and Crane proposed their Riemannian model in p], two specific 
q-( choices of normalization factors for the amplitude were given in ^1] and ^Zj ; we cai l these the DFKR 

normalization and the Perez-Rovelli normalization, respectively. The DFKR normalization is believed 
(3JT), to lead to divergent amplitudes -4(A) for most triangulations A; this has been tested numerically ^ 

and is thought to reflect some residual gauge symmetry |15| . The Perez-Rovelli normalization has 
been proven to give a finite amplitude -4(A) for each triangulation, but it is not known how to 
regulate the sum over all triangulations. For both normalizations, it has been shown that there is 
no interference in the path integral. 

In order to elaborate on the last point, we recall some background. Barrett-Crane spin foam models 
are generally expressed as a sum over certain triangulations A of an amplitude A( A) , and the amplitude 
-4(A) is itself a sum (or integral) over labellings of faces of A by representations of a group. In the 
sum, one restricts to triangulations of 4-manifolds which interpolate between chosen initial and final 
3-geometries. In the case of the Riemannian model, the group is Spin(4) = SU(2) x SU(2), which is 
the double-cover of SO (4), and we restrict to the balanced representations j £§> j which are indexed by 
spins j. Thus the amplitude can be expressed as 

A 

where 

(1.1) A(A) = J^A(A,j f ). 

is 

Here j/ is a labelling of the triangles of A with spins, and the details of A(A,jf) depend on the specific 
normalization chosen. 

What was shown in is that when the quantity A(A,jf) is nonzero, its sign depends only on 
the initial and final geometries (with their spin labellings), and not on the interpolating manifold, its 
triangulation A, or its spin labelling jf. This is the very strong sense in which there is no interference 
in the Riemannian Barrett-Crane path integrals. The conceptual meaning of this result is not fully 



On 
O 

in 
o 

"o 
cr 

> 

X 



1 



POSITIVITY IN LORENTZIAN BARRETT-CRANE MODELS OF QUANTUM GRAVITY 



2 



understood, especially since the model is based on a real-time (e lS ) path integral rather than an 
imaginary-time (e~ s ) path integral. Computationally, the result is extremely powerful, allowing one 
to use statistical methods such as the Metropolis algorithm for calculating expectation values [3]. 

Now we describe the current state of knowledge of the Lorentzian Barrett-Crane model. This model 
was defined by Barrett and Crane in and involves representations of the connected Lorentz group 
SOo(3, 1) from the principal series. These representations are indexed by non-negative real numbers 
p, and so the sum in equation becomes an integral over a product of copies of R + . In |18j . 

a specific normalization was given, which we again call the Perez-Rovelli normalization. With this 
normalization, the Lorentzian Barrett-Crane model was proved ^JEj to give a finite amplitude -4(A) 
for each triangulation. In |H] this proof was extended to include certain degenerate triangulations. As 
in the Riemannian case, it is not known how to regulate the sum over triangulations. 

In [3] numerical computations were done which give evidence to the conjecture that there is no 
interference in the Lorentzian path integral, i.e., that when the quantities A(A,pf) are nonzero, their 
sign depends only on the initial and final geometries. In fact, these amplitudes seem to always be 
non-negative. However, this is still an open question. 

In the present paper, we focus on the case of closed 4- manifolds and show that after transforming 
to a "dual variables" formulation, the model then has the property that there is no interference in the 
new path integral. We also have partial results in the case of manifolds with boundary. 

The dual variables formulation, introduced in |2j , is essentially a specific realization (for Lorentzian 
spin foams) of a transformation first proposed by Pfeiffer |19j in the context of Riemannian spin foams. 
In the dual variables formulation, -4(A) is expressed as an integral with respect to hyperboloid variables 
as follows. For each tetrahedron e of A, there is an associated set of three variables x\ = {x° e , x\, x^). 
Each of these variables takes values on the future 3-hyperboloid. The first comes from the edge 
amplitude, and the other two are associated to the two 4-simplices that e is contained in, and come 
from the vertex amplitude. 

The transformation to dual variables allows us to define an amplitude A(A,x t e ) such that -4(A) is 
the integral of A(A,x l e ) with respect to the x\ variables. In this paper, we prove that the quantities 
A(A,x l e ) are always non-negative. While this doesn't tell us anything about interference between the 
amplitudes A(A,pf), it does imply that in the original path integral there is no interference between 
different triangulations. That is, we find that the quantities -4(A) are always non-negative. 

As we have explained, our result is not exactly analogous to the result in the Riemannian case, in 
that we have not shown that there is no interference between the amplitudes A(A,pf) (although we 
still believe this to be true) but only that there is no interference between the amplitudes A(A,x l e ). 
At first, this would seem to not be enough to allow statistical methods to be used for the computation 
of expectation values of observables, especially observables which are naturally functions of the pf 
variables. But it turns out that in many cases the computation of such an expectation value can also 
be transformed into the dual variables formulation, where there is no interference, and that statistical 
methods can be successfully applied. Using this, the first author has made the first ever computations 
of expectation values of geometric observables in a Lorentzian Barrett-Crane model. These will be the 
subject of a forthcoming paper |1U| . 

In addition to proving the above results for the Perez-Rovelli model, we illustrate the robustness of 
the methods by handling models with the non-standard face normalizations Af = 1/pf and Af — 
The methods could be used to check for positivity in other normalizations as well. 

Next we outline the paper. In Section|21we recall the Lorentzian Barrett-Crane model, describe the 
transformation to dual variables, and show how the results of Section pimply our positivity results. In 
Section|3|we give some elementary results about Fourier transforms that we will need. In Section^Jwe 
prove several positivity results, with each proof corresponding to a different normalization. We give 
conclusions in Section and prove some technical results about Fourier transforms in Appendices lAl 
and|B] 

2. Positivity in the Dual Variables Form of the Perez-Rovelli Model 

In this section, we review the Perez-Rovelli model and its dual formulation in terms of hyperboloid 
variables J3 . Then we show how the results of Section0]can be used to show that there is no interference 
in the dual path integral. As we shall see, in the case of a manifold with boundary, the conclusion 
depends on whether the boundary data is in terms of hyperboloid labellings or representation labellings. 
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For simplicity, we assume throughout that our triangulations are non-degenerate, but this can be 
relaxed in many cases. 

2.1. The Perez-Rovelli Model. We start by recalling the general form of the Lorentzian Barrett- 
Crane model, which assigns to each triangulation A of a 4-manifold an amplitude. For now we deal 
with closed 4-manifolds; we discuss manifolds with boundary in Subsection 12. 31 To be consistent with 
cited papers, we work with the dual 2-skeleton of the triangulation A: to each 4-simplex, we associate 
a dual vertex; to each tetrahedron, we associate a dual edge; and to each triangle, we associate a dual 
(polygonal) face. V, E and F denote the sets of dual vertices, dual edges and dual faces, respectively. 
(This use of the word "dual" has nothing to do with the phrase "dual variables".) 
The general Lorentzian Barrett-Crane amplitude is 



(2.1) 



(2. .2±s\feF J \eeE J Wv / feF 



„2 



where the factors of the form pj arise from the measure on the principal series of SOo(3, 1) representa- 
tions and the amplitudes Af, A e , and A v will be defined shortly. It should be noted that while some 
authors absorb the pj factors from the measure into the definition of Af , in the present work we keep 
the measure and the face amplitude distinct. 

To more compactly represent the multiple integrations that take place in l|2.1(l and in the amplitudes, 
we introduce the symbol J7 to indicate a formal product of symbols such as integral signs or measures. 
When a relation such as / 9 v appears below a product symbol (multiplicative or formal) , the product 
is taken over all of the objects on the left hand side of relation that satisfy the relation; for example, 
ri/g-u denotes a product over all dual faces / such that v is a member of /. 

We now turn to the definition of the amplitudes A v , A e , and Af that appear in l|2.1|) . First we 
need the following notation. We denote hyperbolic space by H 3 = {x £ M 4 | x ■ x = 1 and x > 0}, 
where x ■ y is the Minkowski inner product x ■ y = x a y a — xiyi — x 2 y 2 — ^32/3- If x and y are points in 
hyperbolic space, then 4>{x, y) denotes the hyperbolic distance cosh -1 (a; ■ y) between x and y. If v is a 
dual vertex which is contained in the dual face /, then <j>l is <p(x ei , x e2 ), where e% and e 2 are the two 
edges of the polygon / which meet at v; equivalently, in the triangulation picture, e\ and e 2 are the 
two tetrahedra in the 4-simplex v which share a common triangle /. 

Up to a regularization that will be defined shortly, the Lorentzian Barrett-Crane vertex amplitude 
of |S] is defined as 

(2.2) MPf) = 

where the kernel function K p , is 




(2.3) K P M) = 



sin(p f (f>l) 
Pf smh(0i ) 

As there are ten triangles in the 4-simplex to which a vertex is dual, the integrand is a product of ten 
such kernels. 

Although the expression for the vertex amplitude given above is generally infinite, we adopt the 
usual regularization by fixing the value of one H 3 variable and dropping the corresponding integral 
over H 3 — the answer is independent of the choices 0- 

When dealing with several of these integrals at once, we distinguish the dummy variables of inte- 
gration x e contributed by different dual vertices v by denoting them x v e . 

While the original model of |S] specifies the vertex amplitude A v as given in l|2.2|) . it leaves un- 
specified the edge and face amplitudes A e and Af. The most common choices of these are due to 
Perez and Rovelli. In this model ^Hjj Perez and Rovelli specify Af = 1 and A e = 84(^1,^2,^3,^4), 
where p±, .. . ,pi are the representation variables labelling the four faces containing e. The function 
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04 (pi, P2iP3-,P4) is known as the eye diagram and can be defined as follows: 



(2.4) 



04(Pl,P2,P3,P4) = / K pi{ < P{x,y))K P2 (4,{x,y))K p;i (4,{x 1 y))K Pi {(l){x,y))dx 

_ J_ /" sin(pi0(x, y)) sm(p 2 (f>(x, y)) sin(p 3 0(x, y)) sin(p 4 <ft(s, y)) 

27r 2 7 H 3 pip 2 p 3 p 4 sinh 4 (<?!>(:E,z/)) 

2 f 00 sin(pir) sinfer) sin(p3r) sinfj^r) 
' ~ — dr. 



dx 



lo pip 2 P3P4sinh (r) 

The value is independent of the choice of y € H 3 . For simplicity, we take y to be the origin O of 
hyperbolic space, so r = </>(x, O) is the radial coordinate of x. When dealing with several of these 
integrals at once, we distinguish the dummy variables of integration x and r by denoting them x® 
and r e . 

Inserting iTT^ . J23J), Af = 1 and it2~4"|l into (|2~T|) gives 

3in(p/r e ) \\ 



(2.5) ^(A)= n/ %)(n^)fn/ *^ 1U 

V/eVo /V/Vf VVeei^ Trsmh 2 ^)^ P/ 

fnf n / ^)(n^))) 

As called for by the regularization, for each dual vertex v, a dual edge e$ is chosen and the integral 
with respect to x v e v is omitted. 

We shall henceforth refer to this choice of amplitudes as the Perez-Rovelli Model. 



2.2. Face Factoring and Positivity for Closed Manifolds. Equation 12.5(1 defines the Barrett- 
Crane amplitude as the result of first integrating with respect to the hyperboloid variables x° e (or, 
cquivalcntly, r e ) and x". and then integrating the result with respect to the representation variables 
Pf. This order of integration has provided the context for most theoretical and numerical work to 
date. In HJ, the integrals are shown to be absolutely convergent for all non-degenerate 4-manifold 
triangulations, and so it follows that the order of integration can be reversed. That is, we can integrate 
with respect to the representation variables p/ first, and regard the integration with respect to the 
variables r e and x v e as the path integral. The result is 

(2.6) .4(A) = (tT r dr A (m TT / d < ) M( A '< 

where 




(2J) ^'{UTssmM^^pmi. *><».*.•.>* 




and 

(2.8) F f (p f ,4>l,r e 



(Pf ' ' ' sil % ^(/,dcg v (/))) S ' m (Pf r *(/.l)) • ■ • sin (?>/ ^(/,deg E (/))) 



deg v (/)+deg E (/)- 
P f 



We have factored the integrand into face factors Ff{pf,(j>l,r e ) that each depend on only a single p/ 
variable (but depend on several x v e and r e — (f)(x°,0) variables). In l|2.8|l . deg E (f) denotes the number 
of edges contained in the face / and e(f,i) selects the ith edge contained in the face /. Similarly, 1 
deg v (f) denotes the number of vertices in the face / and v(f, i) selects the ith edge contained in the 
face /. The number of factors in the denominator is two less than in the numerator because of the pj 
in the measure. The integrals J^ c 'Ff(pf, <pf, r e ) dpf, which we shall refer to as integrated face factors, 
can be found exactly in closed form j^j, which is particularly useful for numerical applications. 

As we shall see in Subsection l4.ll an analysis of the face factors in the Fourier domain shows that the 
integrated face factors J Q F/(p/ ,<f>l, r e ) dpf are non-negative when deg v (f) + deg E (f) > 3, which is 

^For closed 4-manifolds, dual faces are closed polygons so deg B (/) = deg v (f). If the 4-manifold has a boundary, the 
dual object to a triangle in the boundary is an "open polygon" with free edges. In this case, deg B (/) = dcg^(/) + 1. 
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always the case for non-degenerate triangulations. It follows immediately that each amplitude A(A, x\) 
is non-negative, and therefore that the amplitude ./4(A) for any triangulation of a closed 4-manifold 
is non-negative. Physically, this means that there is no interference between different triangulations 
A. This result has not previously been shown, although there were strong indications from numerical 
work 4: as pointed out in the introduction. 

The non-negativity of A(A, x\) means that for a given triangulation A, there is no interference 
amongst any of the x\ configurations that can be assigned to its edges. This result is a necessary con- 
dition for applying many statistical mechanical methods such as the Metropolis algorithm to computing 
the expectation values of observables. Remarkably, this result allows one in certain cases to extract 
expectation values of observables depending on the original representation variables (and possibly on 
the hyperboloid variables as well) |l(Jj . 

We note that the face factoring transformation to hyperboloid variables and the positivity results 
are valid for models with choices of edge and face amplitudes different from the Perez-Rovelli model 
(see Subsections 14.21 and I4.3[l . For example, a simplified model with the edge amplitude trivialized 
to A e = 1 and face amplitude taken to be Af = 1/pj is being used to test numerical applications of 
positivity |10| . 

2.3. Positivity for Manifolds with Boundary. We now turn to triangulations of 4-manifolds with 
boundary; in spin foam quantum gravity these are interpreted as histories that interpolate between 
incoming and outgoing boundary 3-geometries. 

The amplitudes for histories are generally required to satisfy a composition law. That is, if Ai is 
a triangulation with initial boundary Ti and final boundary T2, and A2 is a triangulation with initial 
boundary T2 and final boundary T3, then the amplitude for the combined triangulation Ai Yir 2 ^ 2 1S 
required to be the integral of A(Ai)A(A2) over all boundary data on T2- (The dependence on the 
boundary data has been suppressed from the notation.) 

In the original formulation in terms of representation variables, triangles in the boundary tetrahedra 
arc assigned fixed representation labels as boundary data. If this convention is still used, then one 
obtains an amplitude -A(A, {p/}/e&F ) that depends on the labels on the faces in the boundary. The 
composition law is satisfied as long as in Equation l|2.1|l the edge amplitudes A e are replaced with 
\J Ae., when e is in the boundary. (If the face amplitudes are non-trivial, they also must be adjusted in 
the same way.) 

Because of the factors y/Ae, this model cannot be converted completely to face-factored form, as 
the integrals with respect to r e cannot be brought outside of the square root. However, it can be 
partially converted to face-factored form, with some integrals with respect to r e variables left in place. 
An additional difference between this form and Equations (|2.6|l and i|2.7[l is that the integrations with 
respect to pj are omitted for faces / in the boundary. Since the unintegrated face factors l|2.8fl are 
oscillatory functions, there is clearly interference between different x\ configurations in this case. It 
is likely that after integrating out the x\ dependence there is no interference between the resulting 
amplitudes -4(A, {p/} fedF), but our methods do not give any information about this. 

However, since in most cases the majority of the variables lie in the interior of the triangulation 
and hence give non- negative contributions, we are hopeful that statistical methods will be able to take 
advantage of this. In any case, the fact that the integrated face-factors can be analytically computed 
will be very useful for computations. 

2.3.1. A dual variables approach. In the dual variables formulation of the Barrett-Crane model, one 
regards the histories as being labelled by points in hyperbolic space. From this point of view, it is really 
more natural to have boundary data labelled by points in hyperbolic space as well, and to integrate 
over the representation variables pf even for boundary faces. With A e = Q4, the natural boundary 
variables to hold fixed in order to satisfy the composition law are the r e variables. (It does not make 
sense to hold any of the x v e variables fixed, since they are not shared when two triangulations are joined 
along a common boundary.) However, the \J A e factors lead to complex- valued face- factors which we 
are unable to analyze. 

If we simplify the model by setting A e = 1, then there are no boundary variables. The model can 
be converted fully to face-factored form, and our methods show that the amplitudes are non-negative. 
Incidentally, this implies that in the model with A e — 1 and representation labels on the boundary, 
the integral of A(A, {j5/}/gof) over all boundary data is non-negative. Note that for this model to be 
finite, we need Af — 1/P/- 
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3. Fourier transforms 
For a > 0, let Xa be the rectangle function of the interval [—a, a], defined by 

!1 for k 6 (—a, a) 
1/2 for fc = ±o 
otherwise 

and let sinc(t) = sin(t)/t. 

Define the Fourier transform T{f) of a function / 6 by 

/OO 
/(t)e-^. 
oo 

This integral converges to a finite value for each k, and defines a continuous function of k |21l Thm. 
7.5], which is in L 2 (R). 

If / is in L 2 (M), then we define the Fourier transform of / in a more indirect way. The product Xaf 
(essentially, / restricted to the interval [—a, a]) is in and so we can define J^iXaf) as above. 

Since Xaf is a l so m L 2 (R), it follows that Fixaf) is in L 2 (R), and we define the Fourier transform of 
/ to be 

Hf) = lim HXaf), 

a — >oo 

with the limit taken in L 2 (M.). So T(f) is only defined up to equivalence in L 2 (R), i.e., only up to 
changes on a set of measure zero. 

The inverse Fourier transform of an L 1 function / is 

(3.1) ^ X (/)(fe) = f- / f(t)e ikt dt = -F(f)(-k). 

The definition is extended to L 2 functions by the same limiting procedure use to define the forward 
Fourier transform. 

The Fourier transform is an isomorphism from L 2 (M) to i 2 (R), and the inverse Fourier transform 
is its inverse Thm. 7.9]. 

Example 3.1. The rectangle function Xa is in L (R), and it is elementary to check that its Fourier 
transform is given by 

!F(xa)(k) = 2asinc(ai). 

The function asinc(ai) is in L 2 (R) but not L 1 (M). To compute its Fourier transform, we use \3.1\l to 
find 

F~ 1 {TTXa)(k) = —F{Xa)(-k) = asinc(ai). 

It follows that the Fourier transform of a sinc(ai) is given by 

J r (asinc(at))(k) — TTXa{k)- 

We have shown this is true as functions in L 2 (J&), i.e., almost everywhere. One can in fact show that, 
suitably interpreted, this is true for all k. 

We recall the following fact relating the Fourier transform and the convolution. See for example 
[HI Thm. 7.2] and 1 1 p. 478]. 

Lemma 3.2. // 

(1) / and g are in L 1 (]R), or 

(2) / and g are tempered distributions and at least one has compact support, 
then F{f*g)=F{f)F{g). 

We will make use of the following Corollary. 

Corollary 3.3. If F and G are tempered distributions whose Fourier transforms satisfy (1) or (2) 
above, then 

T{FG) = i^HF)*HG). 
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A tempered distribution is a distribution that is bounded by a polynomial at infinity. Any LP 
function, 1 < p < oo, defines a tempered distribution. 

Proof. Let / = T(F) and g = T{G). Then F(f*g) = T{f)T{g) by the lemma. By JO]) . J 7 ~ 1 (f*g) = 
2nJ 7 ~ 1 (f)J : ~ 1 (g) and so J r ~ 1 (J r (F)*J 7 (Gj) = IttFG. Applying the Fourier transform gives the desired 
result. □ 

4. Face Factor Positivity 

In the following subsections, we show that the integrated face factors are non-negative for various 
choices of face amplitude. In this section, we introduce new variables ai and t — the ai are non-negative 
real numbers that correspond to the hyperbolic distances <j>l and the non-negative radial variables r e . 
The pf variable at a face is represented in our proofs by t, symbolizing the "time" domain or original 
domain to which Fourier transforms are applied. 

4.1. Positivity for Ap = 1. This choice of amplitude, which is that of the original Perez- Rovelli 
model, leads to face factors of the form 

n 

F%(ai, . . . , a n , t) = t 2 a t sinc(a 4 t). 

i=l 

We wish to show that the following quantity is non-negative, for all choices of the at > 0: 
7^(oi,...,o n ) = / F 2 (a 1 ,...,a n ,t)dt = - F^(a lt . . . ,a n ,t) dt. 

JO 1 J~oo 

Here we have used that F 2 is an even function of t. 

When n — 3, it is straightforward to perform the integration and show that the result is non-negative. 

We now give a proof that I 2 is non- negative for n > 4. Ignoring the positive factor 1/2, I 2 is equal 
to the value of the Fourier transform 

T ^t 2 Y\ a i sinc(aji)^ (fc) 

at k = 0. 

Our method will be to evaluate the Fourier transform for general k in a way that lets us see that 
the value at k = is non-negative. During the calculation we will take the Fourier transform of L 2 
functions, so we must be careful about evaluation at a specific point (like k = 0). Since n > 4, the 
function t 2 Y[ ai sinc(ait) is in L 1 (R), and so its Fourier transform is continuous in k |^ Thm. 7.5]. 
Thus if we ensure that our indirect computation of the Fourier transform leads to a function continuous 
in k, then the value of that function at k = will be equal to I 2 . 

Our first step is to use Corollary 13.31 as well as the fact that the Fourier transform of tg(t) is 
i-^T{g){k), when T{g) is continuously differentiable. We obtain 



T\t 2 



Y[ at sinc(a.^) J (k) = i-^^\tJJ_ sxnc(ait)J (k) 



t{ J| a; sinc(M) j (k) 

7T d 2 n 



dk 



Xa,(k). 



2"- 1 dk 2 i-. 

Corollary 13.31 applies since the Xa;' s are L 1 functions. Also note that since n > 4, * Xai(k) has a 

i— 1 1 

continuous second derivative. This justifies the first two equalities, and means that evaluating the last 
expression at fc = makes sense. 

By Appendix ^ * Xa { is an even bump, so it has a maximum at the origin. Hence we have 

i— 1 1 

.if ' 
dk 2 i i 

So we have shown that 2^ is non-negative. 



, . * X ai > o. 
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4.2. Positivity for Af = This choice of amplitude leads to face factors of the form 

a 

Fn(ai,...,a n ,t) = *JJaj sinc(a^)- 

8=1 

We wish to show that the following quantity is non-negative, for all choices of the a, > 0: 

/>OC 

7*(ai, . . . ,a„) = / i^(ai, . . . , a n ,t) dt. 



For n = 2, these integrals converge conditionally when ao =/= a%. In this case, it is straightforward to 
perform the integration and show that the result is non- negative. However, when ao = a%, the integral 
is divergent. Because the n = 2 case does not arise for any non-degenerate triangulation, the ao = a\ 
divergence is not a problem. 

We now proceed to prove positivity for n > 3. To apply Fourier methods, we express as an 
integral over the real line, by multiplying by the Heaviside function 




So we wish to show that 

/oo n 
H(t)t JJflj smc(ait)dt 
00 i=i 

is non-negative. This is the value of the Fourier transform 

T ^H(t)tf[a l smc(a l t) S j (k) 

at k = 0. 

As in the Subsection l4.ll our method will be to evaluate the Fourier transform for general k in a way 
that lets us see that the value at k = is non-negative. Since n > 3, the function H (t) t Y[ di sinc(aii) 
is in Z-^R), and so its Fourier transform is continuous in k. We will ensure that our computation of 
the Fourier transform leads to a function continuous in k, so that the value of that function at k = 
will be equal to 

Proceeding as in Subsection 14. II we obtain 



^"(-ff(t)t jjaisinc(ait) j(k) = —J r (H)(k)*J r (t^a i smc(a i t)\(k) 
( 4 -2) = ^(H)(k)*i^r(^(f[a i sm C (a i t))(k) 



27T v /w dfc 

-^(^)-p^)*4(i!i x - (fc) 

Some comments are in order. First, since n > 3, Lemma IA.5I shows that * x ai is continuously 

i=l 1 

diffcrentiable, which justifies the second equality. Second, the Heaviside function H is not in L P (R) 
for any p, and so we must regard it as a tempered distribution. Its Fourier transform involves S(k) 
and pvi. Here 5(k) is the usual delta distribution and pv^ is the "principal value of jr" distribution 
(see Appendix 0. Third, since F{H) is tempered and J-(t J\7=i ai sinc(aji)) has compact support, 
Corollary 13. 31 justifies the first equality in Q4.2JI . 

Observe that the integral (|4.1|l is real, and so we can ignore the imaginary part of l|4.2|l . (One can 
also use the methods below to see directly that the imaginary part is continuous in k and vanishes at 
k = 0.) So we focus on the real part of the right-hand side. Ignoring positive constants, this is 

(4.3) w\*9'{k) where g(k) = .* Ya^k). 

rC i — 1 



POSITIVITY IN LORENTZIAN BARRETT-CRANE MODELS OF QUANTUM GRAVITY 9 

By Corollary IA.2I g is an even bump. Moreover, by Lemma fA. 51 g' exists and is itself piecewise C 1 . 
So by the results of Appendix FBI the convolution pv^ * g'(k), which is a priori a distribution, is in fact 

the continuous function pv r_°° 9 d£. Thus we can sensibly evaluate this at k = 0, which gives 




But since g is an even bump, g'(—£) is non-negative when £ > and non-positive when £ < 0. Thus 
the integrand is non-negative for all £ ^ 0, and so the result is non-negative. 

4.3. Positivity for Ap = This choice of amplitude leads to face factors of the particularly simple 
form 

n 

F°(ai, ...,a n ,t)= JJa,sinc(ait). 

i=l 

The integrated face factors are 

/>oo /■OO n 

I°(ai, ■ ■ ■ , a n ) ee / F°(ai, . . . ,a n ,t) dt = / JJ a; sinc(a.^) dt. 
Jo Jo i=x 

These are shown to be non-negative, for all n > 1 and all a* > 0, in 0, using methods similar to those 
used here. In addition, they give a closed form for /°(ai, . . . , a n ) and relate it to the volume of certain 
polyhedra and their intersections. See also |H] for an introduction to this material. 

5. Conclusions 

We have shown that in the path integral for the Lorentzian Barrett-Crane model, with various 
normalizations, the amplitude -4(A) for a triangulation A of a closed 4-manifold is always non-negative. 
This provides further evidence for the conjecture 0] that the amplitude A(A,pj) for a triangulation 
labelled by representations is always non-negative, as is the case for the Riemannian model. Our 
method is to use Cherrington's face factoring approach to express the amplitude in terms of hyperboloid 
variables x\ instead of representation variables, and to show that A(A,x l e ) is non-negative. 

The above results hold for the original Perez-Rovelli model as well as those with non-standard face 
normalizations Af — l/pf and Af — l/pj- The methods could be used for other models as well. 

The non-negativity of the amplitudes A(A,xl) is highly significant from the point of view of nu- 
merical computations, as it allows statistical methods such as the Metropolis algorithm to be used for 
efficient computation of expectation values of observables. Moreover, exact expressions for A(A,x l e ) 
can be found |§] . Computational results will be presented in a forthcoming paper _10_ . 

In the case of a 4-manifold with boundary, we find that there is no interference in the interior of 
the spin foam, which may still be useful from a computational point of view. 

Acknowledgements. The authors would like to thank Igor Khavkine, Gord Sinnamon and Josh 
Willis for helpful conversations. 

Appendix A. Convolutions of Xa's 

Since the Fourier transform of sinc(ai) involves the rectangle function \a of the interval [—a, a], and 
we want to understand integrals of products of sine functions, we will study convolutions of the \ a 
functions. A key property of Xa that we will use is that it is an "even bump" . 

Definition. An even bump is a function e : K — > M such that 

(1) e(x) — e(—x) for all x (e is an even function). 

(2) e(x) > for all x. 

(3) e has compact support. 

(4) e(x) > e(y) for < x < y. 

It follows that e is bounded above by e(0) and is in L 1 (R). 



Lemma A.l. The convolution of an even bump with an even bump is again an even bump. 
We thank Gord Sinnamon for the elegant proof of condition Q). 
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Proof. Let e and / be even bumps, and define g to be the convolution: 

/>OC 

g(x) = (e*f)(x) = / e{y)f{x-y)dy. 



This makes sense for every x, since e and / are bounded and have compact support. 
We shall establish each of the four properties for g in turn. 

(1) Using that e and / are even, and then substituting y = —y', we have 



oc 



g(x) = / e{y)f{x- y)dy = I e{-y)f{-x + y)dy=l e{y')f{~x~y')dy'=g{-x). 



(2) Since e and / are non-negative, their convolution will clearly be non-negative. 

(3) Since e and / have compact support, their convolution will also have compact support. 

(4) It suffices to show that for x > and s > 0, g(x — s) > g(x + a). Evaluating the difference, we 
obtain 

g(x — s) — g(x + a) = (e * f)(x — s) — (e * f)(x + s) 
e (y)(f( x - s-y)- f(x + s-y)) dy 

e(x - y')(f{y' - s) - f(y' + s)) dy' (y' = x - y) 

D 

e(x - y')(f(y' - a) - f(y' + s)) dy' + / e(x + y')(f(-y' - s) - f(-y' + s)) dy' 



(e(x - y') - e(x + y'))(f(y' - s) - f(y' + a)) dy' . 

In the last integral, x, y' and s are all non-negative. Thus \x— y'\ < \x+y'\ and \y' — s\ < \y' + s\, 
so it follows from properties (1) and (4) that e(x — y') > e(x + y') and f(y' — s) > f{y' + a). 
Therefore the last integrand is a non-negative function. So g is a non-increasing function on 
the positive axis. 

We have shown that g is an even bump. □ 

Corollary A. 2. An n-fold convolution \a t * • • • * Xa„ of rectangle functions Xa t is an even bump. 

Proof. Since a single rectangle function \a is an even bump, this follows inductively from Lemma lA.ll 

□ 

We also need to make use of smoothness properties of the n-fold convolution of x a 's. 

Lemma A. 3. Let f G L 1 (R). Then f * Xa is a continuous function. If f is continuous, then f * Xa 
is differentiable, and 4j:(f * Xa)(k) is the continuous function f(k + a) — f(k — a). 

If follows that if / is C™ (i.e., / is n-times continuously differentiable) then f * Xa is C" +1 . 

Proof. Suppose / is in i 1 (R). Then 

poo pk+a 

(f*Xa){k)= f(£) X a(k-£)de= f(£)d£, 



which is continuous as a function of k. If / is continuous, then 



d ,. d 



k-\-a 



rffc -(/*Xa)(fc) = ^ / f(£)de = f(k + a)-f(k-a). 



□ 



Definition A. 4. A function / : M — > R is piecewise C 1 if the real line can be partitioned into closed 
intervals [s,t] such that f'(k) exists on (s,t), the one-sided derivatives exist at the end-points, and 
all of these taken together form a continuous function on the closed interval [s,t]. It follows that / is 
continuous on R. 
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Lemma A. 5. Consider the n-fold convolution g — * Xa»- If n = 2, then g is piecewise linear. If 

i— 1 1 

n > 3, then g is (n — 2) -times continuously differentiable, and the (n — 2)nd derivative is piecewise 
linear. In particular, for n > 3, g' is piecewise C 1 since it is either C 1 or piecewise linear. 

Proof. If n = 2, one can check directly that g is piecewise linear. The second part follows by induction 
from the first part and the previous lemma. □ 

Appendix B. Principal value of 1/k 

The Fourier transform of the Heaviside function H involves the "principal value of 1/ fc" distribution 
denoted pvjr. This distribution is defined by 

< P .l,, sp ,.£f). s! »(£ + f)f),, 

This makes sense when <f> has a continuous first derivative and compact support [111 Example 1, page 
439]. 

For the following, refer to the definition of piecewise C 1 from Appendix lAl 

Theorem B.l. If f is piecewise C 1 and has compact support, then the distribution pv^ * f(k) comes 
from the function pv J^°° ^^T^ dk, and this function is continuous. 

In our applications, / is the derivative of a convolution of at least 3 rectangle functions, so by 
Lemma FA. 51 it is piecewise C 1 . 

Proof. Since pvir is tempered and / has compact support, the convolution is well-defined. It is the 
distribution such that 

(pv~ * f{k),4>) = <pv-j- x f(t), cj>{k + £)) = (pv-j-, {f{£), <t>{k + £))) 

for any test function <j> (smooth, with compact support). Note that (/(€), <j>(k + £)) is also a test 
function as a function of k. Expanding the right-hand side gives 

(pvp r f{i)<t>{k+i)dt) 



(pvp f } fi£-k)ct>{l)< 



-OO 

f OO 



pv / i / fit- k)<j){£)d£dk 

-oo * J— oo 

=hm ( r + n r^- k)m d£dk 



k 



lim / | / +/ ] l^m d kd£ 



e— >0 J — co 

/oo 
<f>{£)F e i£) d£, 
oo 

where F e (£) = ^Jl^ + J^°°J /(^ — k)/kdk. The exchange of integration in the second-last equality 
is justified because the integrand is continuous with compact support, since a strip around k = is 
excluded. 

By Lemma fB. 21 lim^ + F c (^) exists pointwise. By Lemma fB. 31 the functions <f>F t are bounded by 
an L 1 function. So by Lebesgue's dominated convergence theorem, we can exchange the limit and the 
integration in the last displayed equation. Thus 

(pv-U/(fc),0) = f lim F e {t)4>i£) d£ = I P v / ~^ dk<t>{£) d£. 
So the distribution pv-jj * f{k) comes from the function P v J^^^— ' dk = lim^ a +F e {£). 
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It remains to show that this function is continuous. Since each F e is continuous, it suffices to show 
at 
that 



that the convergence to lim. F e is uniform. The method used in the proof of Lemma IB. 31 shows 



|lim FS)-F s (£)\<2Ssup\f'\, 

where sup |/'| is defined there as well. 

We have shown that the distribution pv^*/(fc) comes from the continuous function pv c ^ c ' 

□ 

Lemma B.2. Let f be a function with compact support such that both one-sided derivatives exist 
everywhere. Then, using the notation above, lim^ +F e (£) exists for each £. 

Proof. Fix £ and choose A large enough so that the support of f(£ — k) as a function of k lies in [—^4, A]. 
(A will depend on £.) By symmetry, 

Therefore 



lim FJl) = lim ( [ \ D l^—tldk 



im, II I ) f -^-r^-dk 



-A 



k 



= lim | / + / \W- k )-^ dk . 



The last integrand is clearly continuous as a function of k on (0, A]. Since / has one-sided derivatives 
everywhere, it extends to a continuous function on [0, A]. Thus 



A 



inn / W- k \-W dk 



lim / w-k)-m dk 



exists. Similarly, 



exists. □ 

Lemma B.3. Let f be piecewise C 1 with compact support and let <fi be a continuous function with 
compact support. Then, there is a constant E such that for all e > 0, 4>(£)F e (£) < E(f>(£). In particular, 
the family of functions (j>F e is dominated by an L 1 function. (The notation F e is from the proof of 
Theorem \B.1\ and is recalled below.) 

Proof. Choose B and C so that the support of / is contained in [— B, B) and the support of cf> is 
contained in [-C, C]. Let A = B + C. Then for all e > and all £ in [-C, C], 

dk 



k 

A \ f(l-k) 
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Therefore, 

\Fe(£)\ < 




f(t -k)- f(£) 



dk 



rA> 

-[J A + l, l SUPl/ ' |dfc 

<2Csup|/'|, 

where the supremum is taken over R (or [— C, C]) and |/'| means the maximum of the absolute values 
of the one-sided derivatives of /. This supremum is finite, since / is piecewise C 1 and has compact 
support. 

Since this bounds F e (£) for I in the support of 0, we get that <f>(£)F e (£) < 2Csup \f'\<f>(£) for all 
1. □ 
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